Effective three-body interactions in triangular optical lattices 
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We demonstrate that a triangular optical lattice of two atomic species, bosonic or fermionic, can 
be employed to generate a variety of novel spin- 1/2 Hamiltonians. These include effective three- 
spin interactions resulting from the possibility of atoms tunneling along two different paths. Such 
interactions can be employed to simulate particular one or two dimensional physical systems with 
ground states that possess a rich structure and undergo a variety of quantum phase transitions. 
In addition, tunneling can be activated by employing Raman transitions, thus creating an effective 
Hamiltonian that does not preserve the number of atoms of each species. In the presence of external 
electromagnetic fields, resulting in complex tunneling couplings, we obtain effective Hamiltonians 
that break chiral symmetry. The ground states of these Hamiltonians can be used for the physical 
implementation of geometrical or topological objects. 
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I. INTRODUCTION 

With the development of optical lattice technology 
uL Ll 01 , considerable attention has been focused on the 



experimental simulation of a varietyof many-particle sys- 
tems, such as spin chains [JiLAlIlI- This provides the 
possibility to probe and realise complex quantum models 
with unique properties in the laboratory. Such examples, 
that are of interest in various areas of physics, are the 
systems that include many-body interactions. The latter 
have been hard to study in the past due to the difficulty in 
controlling them externally and isolating them from the 
environment Q . To overcome these problems, techniques 
have been developed in quantum optics j^, llOl UA\ which 
minimise imperfections and impurities in the implemen- 
tation of the desired structures, thus paving the way for 
the consideration of such "higher order" phenomena of 
multi-particle interactions. Their applications could be 
of much interest to cold atom technology as well as to 
condensed matter physics and quantum information. 

In this paper we obtain the interaction terms of bosonic 
or fermionic lattices of two species of atoms, denoted here 
by T and | (see 0, S 0). These can be two different 
hyperfine ground states of the same atom coupled via 
an excited state by a Raman transition. The system 
is brought initially into the Mott insulator phase where 
the number of atoms at each site of the lattice is well 
defined. By restricting to the case of only one atom per 
site it is possible to characterise the system by pseudo- 
spin basis states provided by internal ground states of the 
atom. Interactions between atoms in different sites are 
facilitated by virtual transitions. These are dictated by 
the tunneling coupling J from one site to its neighbours 
and by collisional couplings U that take place when two 
or more atoms are within the same site. 

In the following we consider the case of weak tunnel- 
ing couplings, J <C U , assuring that we are always in the 



Mott insulator regime. Our aim is to construct a pertur- 
bative study of the effective interactions with respect to 
the small parameter J/U. Up to the third order this ex- 
pansion will provide Hamiltonians that include three-spin 
interactions. These multi-particle interactions can be, in 
principle, realised with the near future technology. The 
main physical requirement is large collisional couplings, 
U, which can be obtained experimentall y b y Feshbach 
resonances [H llllll- First theoretical O and exper- 
imental [2(| advances are already promising. Hence, the 
time interval needed for those higher order terms to have 
a significant effect can be well within the coherence times 
of the system. 

Several applications spring out from our studies. The 
systematic description of the low energy Hamiltonian 
provides the means for the advanced control of the three- 
spin interactions simulated in the lattice. Hence, different 
physical models can be realised, with ground states that 
present a rich structure such as multiple degeneracies 
and a variety of quantum phase transitions 0, HIH Il4| . 
Some of these multi-spin interactions have been theo- 
retically studied in thepast in the context of the hard 
rod boson 13 > using self-duality symmetries 

|26| . Phase transitions between the corresponding 
ground states have been analysed [2?], . Subsequently, 
these phases may also be viewed as possible phases of 
the initial system, that is in the Mott insulator, where 
the behaviour of its ground state can be controlled at 
will 0. 

The paper is organised as follows. In Section [H] we 
present the physical system and the conditions required 
to obtain three-body interactions. The effective three- 
spin Hamiltonians for the case of bosonic or fermionic 
species of atoms in a system of three sites on a lattice 
are given m Sec. Hill In Sec. |lV]we study the effect Ra- 
man transitions can have on the tunneling process and 
generalised effective Hamiltonians are presented that do 
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not preserve the number of atoms of each species. These 
are of particular interest for the construction of certain 
geometrical evolutions. In Sec. ^complex tunnelings are 
considered and the generation of chiral ground states is 
presented. In Sec. IVII our results are extended toward 
the construction of one dimensional models and several 
applications are discussed. In Sec. IVIII we present an out- 
look and the conclusions. Finally, in the Appendix, two 
alternative methods are presented for the perturbation 
theory that results in the three-spin interactions. 

II. THE PHYSICAL MODEL 

Let us consider a cloud of ultra cold neutral atoms su- 
perimposed with several optical lattices 0, IE IE IE E3 • 
For sufficiently strong intensities of the laser field this sys- 
tem can be placed in the Mott insulator phase where the 
expectation value of only one particle per lattice site is 
energetically allowed ,3'j . We are interested in the partic- 
ular setup of lattices that form an equilateral triangular 
configuration, as shown in Fig. ^ This allows for the 
simultaneous superposition of the positional wave func- 
tions of the atoms belonging to the three sites. As we 
shall see in the following this results in the generation of 
three-spin interaction. 




FIG. 1: The basic building block for the triangular lattice 
configuration. Three-spin interaction terms appear between 
sites 1, 2 and 3. For example, tunneling between 1 and 3 can 
happen through two different paths, directly and through site 
2. The latter results in an exchange interaction between 1 
and 3 that is influenced by the state of site 2. 



The main contributions to the dynamics of the atoms 
in the lattice sites are given by the collisions of the atoms 
within the same site and the tunneling transitions of the 
atoms between neighbouring sites. In particular, the cou- 
pling of the collisional interaction for atoms in the same 
site are taken to be very large in magnitude, while they 
are supposed to vanish when they are in different sites. 
Due to the low temperature of the system, this term is 
completely characterised by the s-wave scattering length. 
Furthermore, the overlap of the Wannier wave functions 
between adjacent sites determines the tunneling ampli- 
tude, J, of the atoms from one site to its neighbours. 
Here, the relative rate between the tunneling and the 
collisional interaction term is supposed to be very small, 



i.e. J <§; U, so that the state of the system is mainly 
dominated by the collisional interaction. 

The Hamiltonian describing the three lattice sites with 
three atoms of species a = {T,!} subject to the above 
interactions is given by 

H = H ( - a) +V, (2.1) 

with 

H w = i ^2 Uaa' ■ njorijv' :, 

jcrcr' 

^=-^(J>t ff a J+1(7 +H.c.), 

3° 

where a JCT denotes the annihilation operator of atoms of 
species a at site j. The annihilation operator can de- 
scribe fermions or bosons satisfying commutation or an- 
ticommutation relations respectively given by 

[aj<T,a\ tal ]± = S jy 5 aa >, 
[aja^fa'U = [a} a ,a), a ,)± = 0, 

where the ± sign denotes the anticommutator or the com- 
mutator. The operator rij a is the corresponding num- 
ber operator and : ... : denotes normal ordering of the 
product of the creation and annihilation operators. The 
Hamiltonian is the lowest order in the expansion 

with respect to the tunneling interaction. 

Due to the large collisional couplings activated when 
two or more atoms are present within the same site, the 
weak tunneling transitions do not change the average 
number of atoms per site. This is achieved by adiabatic 
elimination of higher population states along the evolu- 
tion leading eventually to an effective Hamiltonian (see 
Appendix) . The latter allows virtual transitions between 
these levels providing eventually non-trivial evolutions. 
As we shall see in the Appendix it is possible to describe 
the low energy evolution of the bosonic or fermionic sys- 
tem up to the third order in the tunneling interaction by 
the effective Hamiltonian 

ffrfS _yW + vWk, ( 2.3) 

^ EL ^ E^E S y ' 

The indices a, (3 refer to states with one atom per site 
while 7, 8 refer to states with two or more atomic popu- 
lations per site, E^ are the eigenvalues of the collisional 
part, while we neglected fast rotating terms effec- 

tive for long time intervals. 

It is instructive to estimate the energy scales involved 
in such a physical system. We would like to have a sig- 
nificant effect of the three-spin interaction within the de- 
coherence times of the experimental system, which we 
can take here to be of the order of several 10 ms. It is 
possible to vary the tunneling interactions from zero to 
some maximum value which we can take here to be of 
the order of J/h ~1 kHz Q- In order to have a signif- 
icant effect from the term J 3 /U 2 produced within the 
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decoherence time one should choose U/H ~10 kHz. This 
can be achieved experimentally by moving close to a Fes- 
hbach resonance 0, 0, 0, 0] , where U can take sig- 
nificantly large values as long as trap losses, attributed 
to three body collisions or production of molecules, re- 
main negligible. With respect to these parameters we 
have (J/U) 2 ~ 10~ 2 , which is within the Mott insulator 
regime, while the next order in perturbation theory is an 
order of magnitude smaller than the one considered here 
and hence negligible. Note however, that new interaction 
terms arise only in 5*' 1 order in perturbation due to the 
triangular geometry of the optical lattice. This places 
the requirements of our proposal for detecting the effect 
of three-spin interactions within the range of the possible 
experimental values of the near future technology. 

Within the regime of single atom occupancy per site 
it is possible to switch to the pseudo-spin basis of states 
of the site j given by | j) = = l,riji = 0) and 

| |) = |rtj| = 0, nji = 1). Hence, the effective Hamil- 
tonian can be given in terms of Pauli matrices acting 
on states expressed in the pseudo-spin basis. The sym- 
metries of the initial Hamiltonian H restrict to a large 
degree the form of the low energy expansion. For exam- 
ple, the conservation of the atom number of each species 
corresponds, in the spin basis, to the conservation of the 
total z-spin. Hence, any rotation on the icy-spin plane 
leaves the Hamiltonian invariant. This fact limits the 
possible spin operators that can contribute to the effec- 
tive low energy interactions. Possible terms of the ef- 
fective Hamiltonian are given by for the one-body 
interaction, {cr|cr| +1 , crjcrj +1 +< T j< J j+i} f° r the two body 
interaction or {crf<7f +1 o-| +2 , OJct* +1 + a v cj y J+1 )a z +2 } for 
the three-body interactions where 04 = ui, (see Fig. 0. 
As we can easily verify, the three-spin operators break 
parity symmetry, which is explicitly given by the trans- 
formation T^-l- This indicates that their coupling co- 
efficient should also be asymmetric with respect to this 
transformation, as the original atomic system possesses 
this symmetry. Indeed, in the next section we shall see 
how these terms are generated in the optical lattice setup. 

Another important insight for the ground states of the 
presented Hamiltonians comes from the geometry of the 
lattice. In the case considered here, the triangular pat- 
tern allows for the generation of exotic ground states due 
to frustration, that is ground states that are not minimis- 
ing the energy of the individual Hamiltonians of each link 
of the triangle. This effectively allows for the presence of 
multiple degeneracy in the ground state of the system as 
we shall see in particular examples. 



the lattice without the application of any external field. 
Different rates in the tunneling parameter can then be 
achieved by tuning the intensities of the laser field cor- 
responding to the different directions of the triangle. By 
applying the perturbative expansion (|2.3|l up to the third 
order we obtain that the system can effectively be de- 
scribed by 
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A( 3 Vjaf +1 <7? +2 +Af (aJaJ +1 a? +2 +aJaj +1 ajf +2 ) 

(3.1) 

where <r" is the a Pauli matrix with the usual commu- 
tation properties [cj,cr£] = 2ie ufJ,u '5jk(Jj' ■ The three-spin 
interactions presented in the last line can be viewed as 
the two spin interactions of the second line controlled by 
the third spin (being spin up or down) through an ad- 
ditional a z operator. The couplings, A, B, and A**' , are 
given as expansions in J a /U aa > by 
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III. THE EFFECTIVE THREE-SPIN 
INTERACTIONS 

A. The bosonic model 

Consider the low energy evolution of the triangular 
system given in Fig. ^ of three atoms in three sites of 



where the symbol (T^j) denotes the repeating of the 
same term as on its left, but with the j and J, indices in- 
terchanged. The A term contributes to an overall phase 
factor in the time evolution of the system and can be 
ignored. The B term can easily be eliminated and an 
arbitrary magnetic field term of the form J^j B ' & can De 
added by applying a Raman transition with the appropri- 
ate laser fields. The behaviour of the effective couplings 
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as functions of the tunneling and collisional couplings is 
given in Fig. [21 



(a) (b) 




FIG. 2: The effective couplings (a) A (1) , (b) A (2) , (c) A (3) 
and (d) A' 4 ' as functions of the tunneling couplings J T /U 
and J^/U, where we have set the tunneling couplings to be 
Ji — J% = J£ and the collisional couplings to be Uff = fTj-j = 
Un — U. All the parameters are normalised with respect to 
U. 



One can isolate different parts of Hamiltonian l|3.1f) . 
each one including a three-spin interaction term, by vary- 
ing the tunneling and/or the collisional couplings appro- 
priately so that particular terms in i|3.1fl vanish, while 
others are freely varied. An example of this can be seen 
in Fig. |3 where the couplings A'- 1 -' and A' 3 ' are depicted. 
There, for the special choice of the collisional terms, 
Uff = Un = 2.12Ufi, the A^ coupling is kept to zero for 
a wide range of the tunneling couplings, while the three- 
spin coupling, A^ 3 -* , can take any arbitrary value. One 
can also suppress the exchange interactions by keeping 
one of the two tunneling couplings zero without affecting 
the freedom in obtaining arbitrary positive or negative 
values for A^ 3 ) as seen in Fig. |3] 



Hence, the one dimensional Hamiltonian of the form 
H(B X ,B Z ) = -J2 (B x ** + B z a] + a?a* j+1 a* j+2 ) 

3 

can be simulated in the optical lattice where all of its cou- 
plings can be arbitrarily and independently varied. The 
three-spin interaction term of this Hamiltonian possesses 
fourfold degeneracy in its ground state, spanned by the 
states {| TTT>, I TU),| IU>,| lit)}- The criticality be- 
haviour of this model has been extensively studied in the 
past 0, , where it is shown to present first and sec- 
ond order phase transitions. In particular, for B z — 




FIG. 3: The effective couplings A (1) and A (3) are plotted 
against J T /U and J 1 /U for U n = U u = 2.12C/ and U u = U. 
The coupling A' 1 ' appears almost constant and zero as the 
unequal collisional terms can create a plateau area for a small 
range of the tunneling couplings, while A' 3 ' can be varied 
freely to positive or negative values. 

its self dual character can be demonstrated [25l |2fj . To 
explicitly show that let us define the dual operators 

oc 

~ X rr Z rr Z ^ Z ^- Z I I „X 

k=0 

that also satisfy the usual Pauli spin algebra. We can 
re-express the Hamiltonian H(B x ,0) with respect to the 
dual operators obtaining finally 

H(B x ,0) = B x H(B~\0). 

This equation of self duality indicates that if there is one 
critical point then it should be at \B X \ = 1 as has also 
been verified numerically. Furthermore, the two spin in- 
teractions ct|ct| +1 has a degeneracy with a Z 2 symmetry 
while <Jj<Jj +1 <Jj + 2 has a three fold degeneracy leading to 
a Z3 symmetry. By varying the corresponding couplings 
of the effective Hamiltonian it is possible to induce tran- 
sitions to and from the Z 2 and Z3 ordered phases in a 
similar fashion as has been theoretically demonstrated in 



B. The fermionic model 

Alternatively, one can consider the case of fermionic 
atoms and derive the effective interactions they induce 
up to the third order. Compared to the couplings in 
the bosonic case we now have U-\i = U being the only 
one that is present. The Pauli exclusion principle can be 
signalled by having U^f, — * 00 that eventually forbids 
two fermionic atoms of the same species from occupying 
the same site. Keeping terms up to the third order in 
Jj /U and employing the anticommutation relations l|2.2|) 



5 



we obtain the effective Hamiltonian 

3 
3 = 1 



quantum optics |j, |7| . An interesting example involves 
the application of Raman transitions during the tunnel- 
ing process. These transitions involve the direct coupling 
of the two atomic states | and |. Consequently they are 
not atom-number preserving for each of the species. 



where the effective couplings are a function of the initial 
variables of the Hamiltonian (|2.1(1 and 



4" 

„(3) 



(4) 



(2) 



2JJ > 3 +J 3 >> 

1 I f f f _ 7-1 ji) 

2JJ2 V 1 2 J 3 J 1 J 2 J 3/' 



J T J X 



3 
21P 



^i^i+i J l+2 JjJj + iJj +2 )- 



The dependence of the coupling terms on the parameters 
of the initial Hamiltonian is simpler than in the bosonic 
case. Nevertheless, they can express a similar behaviour 
as can been seen in Fig. 0J 




J'/U 



FIG. 4: The effective couplings (a) (b) ^ (2) , (c) /i (3) 

and (d) ^ as functions of the tunneling couplings J T /U 
and J^/U, where we have set the tunneling couplings to be 

Jl = J 2 =Ji- 



If the tunneling constants do not depend on the pseudo- 
spin orientation then any three-spin interaction vanishes. 
Nevertheless, when the tunneling amplitudes depend on 
the spin and by having just one of the orientation with 
non-zero tunneling, then just the diagonal two- and 
three-spin interactions remain. 



A. The general case 

Let us first consider the case where the atoms are 
strongly trapped by an optical lattice as in the previous 
sections. If the lasers producing the Raman transition are 
forming standing waves it is possible to activate tunnel- 
ing transitions of atoms that simultaneously experience 
a change in their internal state. As we shall see in the 
following the resulting Hamiltonian is given by an SU (2) 
rotation applied to each Pauli matrix of the Hamiltonian 

In particular, we shall consider the case of activating 
the tunneling with the application of two individual Ra- 
man transitions. These transitions consist of four paired 
laser beams L%, L 2 and L[, L' 2 , each pair having a blue 
detuning A and A', different for the two different tran- 
sitions. The phases and amplitudes of the laser beams 
can be properly tuned so that the first Raman transition 
allows the tunneling of the state 

|+) = cos0|a) +sin0e-^|6) 

with tunneling rate J + between two neighbouring sites, 
while the second one activates the tunneling of the state 



-i4> 



\b), 



by an additional phase difference of it between the lasers 
L[, L' 2 , with an effective tunneling rate J_. In the above 
equations <f> denotes the phase difference between the Li 
laser field, while tan 6* = |02/f2i|, where Vti are their cor- 
responding Rabi frequencies. Hence, the effective tunnel- 
ing term is given by 



Vc 



-t 



-i+i 



t 



H.c. 



where the tunneling couplings J+ and J_ are given by the 
potential barrier of the initial optical lattice superposed 
by the potential reduction due to the Raman transitions 
In addition, the creation and annihilation operators are 
given as an SU(2) rotation of the initial ones, i.e. 



9( 



with the unitary SU (2) matrix 



IV. RAMAN ACTIVATED TUNNELINGS 



ff(<M) 



cos U sin V 
sin 9 -e^cost 



A number of variations of the previous Hamiltoni- 
ans are possible by employing techniques available from 



Hence, the resulting tunneling Hamiltonian can be ob- 
tained from the initial one via an SU (2) rotation V c = 



6 



gVg^ , where the corresponding tunneling couplings are 
formally identified, i.e. J + = ,P and J~ = J*. Note 
that the collisional Hamiltonian is not affected by the 
Raman transitions, and hence it is not transformed un- 
der g rotations. 

It is easy to derive the effective Hamiltonian for this 
transformation using the perturbative expansion. In- 
deed, from expressions (|A3|> and l(A4(l of the Appendix 
we straightforwardly obtain the second and third order 
terms of the Hamiltonian H c g that appear after the ap- 
plication of the Raman transition. They are given by an 
SU(2) rotation that acts on the Pauli matrices of the ini- 
tial effective Hamiltonian. Actually this statement holds 
in all orders of the perturbation theory and reads in its 
generality 



where n is the order of the perturbation. Note that this 
useful result holds not only for the <f> rotations, but also 
for the 9 rotations which, in general, do not commute 
with the collisional Hamiltonian flw . 



00 „t/ 



B. Rotated anisotropic XY model 

We now show that the above presented Raman tran- 
sitions can be employed to obtain, for example, the 
anisotropic XY model. The direction of anisotropy is 
determined by the phase difference of the laser fields em- 
ployed for the Raman transition. In particular consider, 
as in the previous, three sites of the optical lattice in 
a equilateral triangular configuration. For simplicity we 
assume J + = J_ = J and = Un = U^i = U . Then 
the effective Hamiltonian to the third order becomes the 
rotation 



H eS {<l),Tr/2)=g{<j))H eS g\<l>), 



(4.1) 



where g(4>) = g{<fi,9 = 0) is a z-axis rotation and H c h 
is the 9 — 7r/2 rotated effective Hamiltonian around the 
y-axis given by 

H cS = ^ (AI + Bat + + ^ (3) «+i^ +2 ), 



with 



3 J J 3 „ J 2 11J 3 

A = 3^r, B = -2 jr. 

2 U U 2 U 2 U 2 



,0) _ 



2 U 



U 2 ' 



1 J 3 
6 U 2 ' 



These effective couplings agree with the ones presented in 
(|3.2|) . Moreover, by controlling the amplitude of the ini- 
tial standing waves that trap the atoms in their equilib- 
rium positions it is possible to reactivate the tunnelings 



and . This has the effect that the overall Hamilto- 
nian is the sum of the two Hamiltonians, the rotated one 
(14.1(1 and the initial one 1(3. 1[) . 

One can now check that the Hamiltonian 1(3.1(1 is invari- 
ant under g{4>) rotations. On the other hand, when we 
add the Hamiltonians Hq and the one from 1(3.11) we ob- 
tain the generalised version of the anisotropic XY model 
where additional third order terms are present. Hence, 
by turning on the and tunnelings we can obtain the 
rotated version of the anisotropic XY model, where the 
rotation is performed with respect to the z-spin axis by 
an angle <f>. This approach provides a variety of control 
parameters (e.g. the angle </> and the ratio of the cou- 
plings of the two added Hamiltonians) and, in addition, 
one can have these variables independent for each of the 
three directions of the two dimensional optical lattice. 
Particular settings of these structures have been proved 
to generate topological phenomena 0, that support ex- 
otic anyonic excitations, useful for the construction of 
topological memories [3j|. In addition, the possibility of 
varying arbitrarily the control parameters of the above 
Hamiltonians and, consequently, of their ground states 
gives us the natural setup to study such phenomena as 
geometrical phases in lattice systems. Examples will be 
presented elsewhere j^- 



V. COMPLEX TUNNELING AND 
TOPOLOGICAL EFFECTS 

Consider the case where we employ complex tunneling 
couplings [3^| in the transitions described above. This 
can be performed by employing additional characteris- 
tics of the atoms like a charge e, an electric moment jl e 
or a magnetic moment fl m , and external electromagnetic 
fields. As the external fields can break time reversal sym- 
metry, new terms of the form {crJcrj' +1 <T| +2 — ct^ctJ +1 <t| +2 } 
appear in the effective Hamiltonian. In particular, the 
minimal coupling of the atom with the external field can 
be given in general by substituting its momentum by 

p^p + eA(x) + ft m x E(x) + (fl e ■ V)A(x), 

where E is the electric field and A is the vector potential. 
All of these terms satisfy the Gauss gauge if we demand 
that V ■ A — and E{r) cx r/r 3 , hence they can generate 
a possible phase factor for the tunneling couplings. 

The first term results in the well known Aharonov- 
Bohm effect [3(| , while the second one is the origin of the 
Aharonov-Casher effect 37j. The first one requires that 
the atoms involved are charged, which is not possible to 
achieve in the optical lattice setup. On the other hand 
it is plausible to consider the electric or magnetic mo- 
ments of the atoms. Nevertheless the Aharonov-Casher 
effect requires that the magnetic moment of the atom 
moves in the field of a straight homogeneously charged 
line, the latter being technologically difficult to imple- 
ment. Though, resent experiments have been performed 
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that generalise the Aharonov-Casher effect partly over- 
coming the technological obstacles j^. The third case 
involves the cyclic move of an electric moment through a 
gradient of a magnetic field finally contributing the phase 



<t> = / (/x e • V)B ■ ds, 
Js 

to the initial state, where S is the surface enclosed by the 
cyclic path of the electric moment. For example, if j2 e is 
perpendicular to the surface S, taken to lie on the x-y 
plane then a non-zero phase, </>, is produced if there is a 
non- vanishing gradient of the magnetic field along the z 
direction. Alternatively, if f2 e is along the surface plane, 
then a non-zero phase is produced if the z component of 
the magnetic field has a non- vanishing gradient along the 
direction of jl e . Hence, it is possible to generate a phase 
factor contribution to the tunneling couplings J = e*^| J| 
with 



(fi e ■ V)A • dx . 



Here Xi and Xi+\ denote the positions of the lattice sites 
connected by the tunneling coupling J. 

In order to isolate the new terms that appear in the 
case of complex tunneling couplings we should restrict 
to purely imaginary ones, i.e. J r J = ±i|JJ|. Then the 
effective Hamiltonian (12.3(1 becomes 



i+l 



- (2) «^+l+<^ 



i+l/ 



(5.1) 

where ei mn with {l,m, n} — {x,y,z} denotes the total 
antisymmetric tensor in three dimensions and summa- 
tion over the indices l,m,n is implied. The couplings 
appearing in (|5.1|l are given in the bosonic case by 



jT 2 Ji 2 JT 2 + Jl 2 JT 2 jl 2 
A = — h — 1 , B = 2 2—, 

2^u ^TT U±i 



^ - Jt \ jl2 J ' 2 + Jl \ T ( 2 ) = 

T {3) = I— ( + — + (T^i), 



Uti \2U-\i U 



-(4) 



'11 

.J^J 1 / 1 



1 



and in the fermionic case by 



.4 



2U 



(T~i), 



B = r< 3 > = 0, 



r(2) 



JV^ f41 JT 2 ^ - JT 

T = 



By taking U^ — > oo, U^-\ = — Un = I . ./ = — J and 
= J, one can set, in the bosonic case, with the aid of 
Feshbach resonances and compensating Zeeman terms, 
all the couplings to be zero apart from t^K Hence, the 
effective Hamiltonian reduces to 



J2 ff i 



a j X (Tfe , 



(5.2) 



(ijk) 



with a = {a x ,ay,(T z ) and t< 4 ' = \ J\ 3 /U 2 . Remarkably, 
with this physical proposal, the interaction term (|5.2|) can 
be isolated, especially from the Zeeman terms that are 
predominant in equivalent solid state implementations. 
This interaction term is also known in the literature as 
the chirality operator |30|. It breaks time reversal sym- 
metry of the system, a consequence of the externally ap- 
plied field, by effectively splitting the degeneracy of the 
ground state into two orthogonal sectors, namely "+" 
and "— " , related by time reversal, T. These sectors are 
uniquely described by the eigenstates of if e ff at the sites 
of one triangle. The lowest energy sector with eigenen- 
crgy E+ = — 2v / 3t( 4 ) is given by 



l*i%> 



1 



V3 



(I TU) + w| T4T) + ^ 2 l 4TT» 



1/2/ 



V3 



(I Ut> + H [M)+uj 2 \ Hi)) 



(5.3) 



U 



2U 2 



The excited sector, \^ ±1 / 2 )i represents counter propa- 
gation with eigenvalue -E_ = 2v / 3r^ 4 ^ and it is obtained 
from H5.3fl by complex conjugation [3(], OHl, • ^ e would 
like to point out that, to the best of our knowledge, this is 
the first physical proposal where this interaction term can 
be isolated, especially from the Zeeman terms that are 
predominant in equivalent solid state implementations. 
Alternative models employing cold atom technology for 
the generation of topologically non-trivial ground states 
are given in 0, El • 



VI. ONE- AND TWO-DIMENSIONAL MODELS 

It is also possible to employ the three-spin interactions 
that we studied extensively in the previous sections for 
the construction of extended one and two dimensional 
systems. The two dimensional generalisation is rather 
straightforward as the triangular system we considered is 
already defined on the plane. Hence, all the interactions 
considered so far can be generalised for the case of a two 
dimensional lattice where the summation runs through 
all the lattice sites with each site having six neighbours. 

The construction of the one dimensional model is more 
involving. In particular, we now consider a whole chain 
of triangles in a zig-zag one dimensional pattern as shown 
in Fig. In principle this configuration can extend our 
model from the triangle to a chain. Nevertheless, a care- 
ful consideration of the two spin interactions shows that 
terms of the form o~faf +2 appear in the effective Hamil- 
tonian, due to the triangular setting (see Fig. |5J). Such 
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Hamiltonian terms involving nearest and next-to-nearest 
neig hbour interactions are of interest in their own right 
^2, El but will not be address here. It is also possi- 
ble to introduce a longitudinal optical lattice with half 
of the initial wavelength, and an appropriate amplitude 
such that it cancels exactly those interactions generating, 
finally, chains with only neighbouring couplings. 



2 4 6 




13 5 7 



FIG. 5: The one dimensional chain constructed out of equi- 
lateral triangles. Each triangle experiences the three-spin in- 
teractions presented in the previous. 



In a similar fashion it is possible to avoid generation 
of terms of the form afaf +2 + of a\ +2 by deactivating 
the longitudinal tunneling coupling in one of the modes, 
e.g. the | mode, which deactivates the corresponding 
exchange interaction. 

As we are particularly interested in three-spin in- 
teractions we would like to isolate the chain term 
E l (<c r f+iC r f+2 + ^f +1 crf +2 ) from the A^ 4) term of 
Hamiltonian Q3.ip . This term includes, in addition, all 
the possible triangular permutations. To achieve that we 
could deactivate the non-longitudinal tunneling for one of 
the two modes, e.g. the one that traps the | atoms. The 
interaction ofof +1 of +2 is homogeneous, hence it does not 
pose such problem when it is extended to the one dimen- 
sional ladder. With the above procedures we can finally 
obtain a chain Hamiltonian as in where the sum- 

mation runs up to the total number N of the sites. 



VII. CONCLUSIONS 

In this paper we presented a variety of different spin 
interactions that can be generated by a system of ultra- 
cold atoms superposed by optical lattices and initiated 
in the Mott insulator phase. In particular, we have been 
interested in the simulation and study of various three- 
spin interactions conveniently obtained in a lattice with 
equilateral triangular structure. They appear by a per- 
turbation expansion to the third order with respect to the 
tunneling transitions of the atoms when the dominant in- 
teraction is the collisions of atoms within the same site. 
Among the models presented here we specifically consid- 
ered the cr? +1 af +2 interaction, a third order generalisa- 
tion of the rotated inhomogeneous XY model, as well as 
interactions that explicitly break chiral symmetry. These 
models can exhibit degeneracy in their ground states and 



undergo a variety of quantum phase transitions that can 
also be viewed as phases of the initial Mott insulator. 

It is possible to employ quantum simulation techniques 
|3l|. in a similar fashion as for two spin Hamiltonians, 
to generate effective three-spin interactions that are not 
possible to obtain straightforwardly from the optical lat- 
tice system. Hence, a variety of additional Hamilto- 
nians can be obtained by considering manipulations of 
the above three-spin interactions with the application of 
appropriate instantaneous one or two spin transforma- 
tions. The possibility to externally control most of the 
parameters of the effective Hamiltonians at will reenters 
our model as a unique laboratory to study the relation- 
ship among exotic systems such as chiral spin systems, 
fractional quantum Hall systems or systems that exhibit 
high-T c superconductivity |2^, |3fJ. In addition, suit- 
able applications have been presented within the realm 
of quantum computation [3^ | where three-qubit gates can 
be straightforwardly generated from the three-spin inter- 
actions. Furthermore, unique properties related with the 
criticality behaviour of the chain with three-spin inter- 
actions has been analysed in [LH where the two-point 
correlations, used traditionally to describe the criticality 
of a chain, seem to fail to identify long quantum cor- 
relations, suitably expressed by particular entanglement 
measures |42|. 

In conclusion, we have presented a physical model that 
can efficiently simulate a variety of three-spin interac- 
tions. The employed optical lattice techniques give the 
possibility to externally manipulate and control the cou- 
plings of the interactions. The effect of these terms will 
eventually be significant with the improvement of the ex- 
perimental techniques. Importantly, the three-spin inter- 
actions can be isolated from two-spin ones or from pos- 
sible Zeeman terms that are always present in the corre- 
sponding spin systems. This makes the further study of 
their properties an important task for future work. 
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APPENDIX A: PERTURBATION THEORY 

Consider the case of two species of atoms trapped 
in optical potentials forming a triangular configuration 
subject to the Hamiltonian given by l|2.f|) . For simplic- 
ity define the diagonal free Hamiltonian to be given by 
■Hij = EiSij, where Ei is either zero or proportional to 
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Uaa, Utb or U a b- As we have already mentioned we con- 
sider the case where tunneling couplings are much smaller 
than the collisional ones J <C U. Then the evolution of 
the system is dominated by the term H^°\ In fact, when 
we start from a configuration of one atom per lattice 
site, denoted by the subspace M of configurations, and 
activate small tunneling couplings, the change of atom 
number per site would be energetically unfavourable and 
is hence adiabatically eliminated. 

To see this analytically we employ the interaction pic- 
ture with respect to the Hamiltonian obtaining 



Hiijit) = Vij exp 



i(Ei - Ej)t/h 



(Al) 



The evolution operator in the interaction picture is given 
by the time ordered formula 



^/(i,0) = Texp --/ Hi(t')dt' 



H!{t')dt' 



~^fdt' J* dt" J* df'H^Hjit^Hjit'") 



+0{{Jtf) 



(A2) 



Higher orders are negligible as long as times t are con- 
sidered for which Jt remains sufficiently small, while Ut 
is large enough to avoid the accumulation of population 
outside the subspace M. The latter condition is nec- 
essary to exempt fast rotating phase factors appearing 
when performing the above time integrals. These phase 
factors are of the form e tuJt — 1 and 

lim (e lult - l) = lim ( - 2 sin 2 — + isinujt) 



which is zero for u> oc Ei — Ej ^= 0. These conditions are 
in agreement with the previous demands that Jt is very 
small while Ut is relatively large. Hence, we can directly 
calculate each term of the expansion (|A2() without having 
to take into account fast rotating terms. 

The effective Hamiltonian H c g that corresponds to this 
evolution can be obtained by the term proportional to 
time t in the expansion of the evolution operator, i.e. 

Wj(t,0) =1- ^H cS t + 0(t 2 ) 
n 

Consider now the second term on the right hand side of 
l|A2| . This term gives no evolution within the subspace 
M of states as the tunneling Hamiltonian term V moves 
you necessarily out of the M configurations. The third 
term gives (see 



(,-H eff ) a p - - —g — ■ 



(A3) 



where a and (3 indicate states in M, 7 indicates states out 
of M and E are the eigenstates of H^ ', where we have 
used E a = Ep = 0. This gives the usual second order 
effective Hamiltonian presented in detail in 0, Q . Con- 
sider now three sites and the effect of the third term in 
Eqn. l|A"2|l . Finally, we obtain the effective Hamiltonian 
with matrix elements 



7(5 



E~E$ 



(A4) 



With the formulae l|A3(l and l|A4(l one can perform the 
perturbation up to the third order and find the desired 
three-spin interactions (|2.3(l . In practise the evaluation 
of the terms that contribute to the three-spin Hamil- 
tonian is quite simple. The states corresponding to 7 
and 8 include sites with two or three atoms of the same 
or of different species. Hence, E y ,Eg oc U aa >. Next 
you need to consider the different evolutions of the form 
a — > 7 — > 5 — > f3, that populations undertake. The 
tunneling couplings J a are determined by each of these 
transitions and an appropriate coefficient is obtained in 
the case of the bosonic generation or annihilation of two 
atoms of the same species in one site. 



APPENDIX B: ADIABATIC ELIMINATION 

As an alternative procedure it is possible to eliminate 
the fast oscillating term without performing the pertur- 
bative expansion. This elimination is related with the 
adiabatic elimination of the states with two or more 
atoms per lattice site that are separated from the states 
with one atom per lattice site (configurations in M) by a 
large energy gap proportional to U aa -i . In fact if we set a 
decomposition of the three site in terms of basis states of 
the form |«i ji, i 2 , J2', is, J3) where 1,2,3 denote the site, 
and ik and jh denote the number of atoms of species | 
and I respectively in site k we can write the general state 
of the three sites as 

i*W) = E c ^( i )i^; j 2,j 2 ; l3 ,j3). 

ik ,3k 



By employing the Schrodinger equation we can obtain 
the time-differential equations of the coefficients c*-* of 
the form 



where H) k: "' 

l k'Jk 



(Cj k k \H\cf'' i ). It is easy to verify that the 
elements of H with indexes (ikjk) corresponding to states 
that do not belong to M include fast rotating phases 
and, hence, they are zero, i.e. for those states cj* = 
0. This provides a set of linear equations of the form 
y\, „•/ Hl klk c] k = that can be solved, in principle, 
explicitly. In our case, Eqn. (|B1|) has overall 56 equations 
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resulting from the Schrodinger equation with 48 reduced terms of small couplings J < [/ obtaining an alternative 
to a linear system of coupled algebraic equations. This verification of the previous perturbative expansion, 
set can be solved by a computer and then expanded in 
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